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(57) Abstract: The invention refers to a method of 

automatic recognition of musical compositions and 
sound signals, which is used for the ideniificaiion 
of musical compositions and sound signals played 
by radio or TV, or performed in public places. Ac- 
cording to this invention, there is a selection of a de- 
sirably large number of musical compositions and 
sound signals, which we want to identify. In every 
one of these signals an original procedure is applied 
leading to the extraction of a set of characteristics 
which will finally represent a model signal. Subse- 
quently, for the implementation of the recognition, 
the unknown musical composition or sound signal 
is received and digitised. To its digitised version 
the same procedure of extracting a set of charac- 
teristics is applied. These are compared with the 
corresponding sets of the model signals and with 
original criteria it is decided if there is a model sig- 
nal that corresponds to the unknown signal under 
consideration. Moreover, it is decided which model 
signal exactly corresponds to the unknown one. 
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Method of Automatic Recognition of Musical Compositions and 

Sound Signals 

This invention refers to a method of automatic recognition of musical 
5 compositions and sound signals and it is used in order to identify musical 
compositions and sound signals transmitted by radio, TV and/or performed in public 
places. 

During the past, efforts for the development of methods for the automatic 
- recognition of musical compositions and sound signals have been made, that led to 

10 the creation of systems performing this task. However, these methods and the 
related systems manifest low percentage of success&l recognition both for the 
musical compositions and for the sound signals of interest. The introduced method 
offers much better percentage of fiiUy automatic recognition, higher or equal than 
ninety eight percent (98%). 

15 According to this invention, there is a selection of a desirably high number of 

musical compositions and sound signals, which we want to identify. For easy 
reference we will refer to these compositions and signals with the term model 
signals. In every one of these signals an original procedure is applied leading to the 
extraction of a set of characteristics which will finally represent each model signal. 

20 Subsequently, for the implementation of the recognition, the unknown musical 
composition or sound signal is received, in which the same procedure of extracting a 
corresponding set of characteristics is applied. These characteristics are compared 
with the corresponding sets of characteristics of the model signals and, by means of 
a number original criteria, it is decided if one (and which one exactly) of the model 

25 signals corresponds to the unknown signal under consideration. This procedure is 
described in figure 1 . 

It is stressed that, officially, there is no reference in international bibliography 
for a similar method or a relative system. In the world market there are very few 
similar systems which offer a percentage of successful recognition less than sixty 

30 percent (60%). 

The invention is described more thoroughly below: 

First, the whole frequency band from 0 to 11025 Hz is divided to sub-bands 
that are almost exponentially distributed. An implementation of such a division 
presented in Table 1 . 

35 According to this implementation, the whole frequency band from 0 to 1 1025 

Hz is divided in 60 sub-bands. 

Subsequently, each model signal is digitised with a random sampling frequency 
preferably greater than or equal to 1 1025 Hz and a window of 8192 or 16384 
or 32768 sample length, slides on the obtained digitised signal. In every such 

40 window, an adaptive Fast Fourier Transform is applied and the Discrete Fourier 
Transform absolute value is obtained. Next, the frequency domain window is divided 
in sections according to the aforementioned frequency sub-bands choice (see Table 
1) and then, in every such section, all the peaks of the absolute value of the Fourier 
transform are spotted and the greater one is obtained. The value of this peak is 

45 called "section representative". Then the L "representatives" with the greater values 
are spotted, where the value of L may vary from 13 to 30, while the most frequently 
used L value is 20. The indicators of the sections corresponding to these 
representatives, sorted in increasing order, form a vector, which constitutes the 
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"representative-vector^' of the window. The above procedure is repeated while the 
window slides on the whole digitised model signal thus creating all the 
representative vectors for the specific model signal. Notice that, while^ the window 
slides on the model signal, the generated representative vectors often remain 

5 unchanged in two successive windows, successive in the sense that they start in two 
positions differing one sample the one from the other. For this reason, in every 
representative vector we assign a number indicating the number of subsequent 
windows in which the specific vector remained unchanged. For that number we will 
use the name "number of repetitions" of the representative vector. For the set of the 

10 generated representative vectors of each model signal we will use the name "the 
model signal set of representatives". The aforementioned procedure is described in 
figure 2. 

For the identification of the unknown sound signal, which from now on will be 
called the "unknown signal", the following procedure is used: 

15 A part of the unknown signal of length varying from eight (8) to sixteen (16) 

seconds is received, digitised and registered, at least temporarily. At the beginning of 
that unknown signal part a window of length Wj^ = 8192 or Wf = 16834 or 
Wf = 32768 samples is obtained; notice that in any case this window will be of the 
same length with the sliding window which was used for the model signals. In this 

20 window a Fast Fourier Transform is applied and the absolute value of it is obtained. 
Afterwards, all the peaks of the absolute value of the Fourier transform are spotted 
and S copies of these peaks are created. For the creation of every copy of the peaks, 
the positions of the peaks are multiplied with a different coefficient /y , i = 0, I, 
S, which is called "window shift coefficient". Thus, S+1 different groups of peaks 

25 are created. For every one of these groups the following procedure is realised: the 
section to which each peak corresponds according to the aforementioned frequency 
sub-bands division, is spotted (see Table 1). For every section to which at least one 
peak corresponds, the greater peak is kept. The value of this peak is called 
"representative of the section of the unknown signal corresponding to the shift 

30 coefficient /, ". 

Next, the L greater value representatives are spotted, where the value of L is 
the same with the one used for the model signals. The indicators of the sections 
corresponding to these representatives, sorted in increasing order, form a vector, 
which constitutes the "first representative vector of the unknown signal 

35 corresponding to the shift coefficient // ". 

Afterwards, the window slides for £ i samples, where the value of ^ j may 
vary from 0,55 * to 1,9 * samples, with most frequently used value the 
^ 1 = 1,4 * . For the new window position and for every shift coefficient 7) , i = 
0, 1,...,S, (S+1) vectors are computed with the way that described above; each such 

40 vector will be called "second representative vector of the unknown signal 
corresponding to the shift coefficient /y ". The above procedure is repeated for M-2 
windows, where each window starts at a sample having a distance of £j samples 
from the start of the previous one, i = 2, 3, M-l, where the value of M may 
fluctuate between 7 and 13 windows, the most usual value being M = 9. In this way 
45 S+ 1 groups of M representative vectors are obtained; for each such group we will 
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20 



employ the name "group of unknown signal representative vectors corresponding to 
the shift coefficient // 

It must be stressed that for a specific application the values, i = 1, 2,. M- 
1, are not necessary equal, but must be kept fixed throughout the whole procedure. 
The exact number (S+1) of the shift coefficients fj varies fi-om 1 to 15, while their 
values are given by the formula: 



\^(L^ySTEP, if i odd 

1, /// = 0 / = 1,2,...,5, 

\-{il2y STEP, if i even 



where STEP is a parameter expressing the shift step, that usually belongs to the 

interval [0.005, 0.01], the more fi-equently used value being 0.0075. The 

10 identification procedure described so far is depicted in figure 3. 

For the realisation of the unknown signal recognition, each group of unknown 
signal representatives is being compared with elements of the set of representatives 
of each model signal separately. To set ideas, each of the S+l groups of M unknown 
signal representatives is compared with groups of M model signal representatives by 
15 means of the method consisting of the following steps: 

El) If the first representative vector of one group of the unknown signal is 
called V| and the first representative vector of the of the model signal is called Ui , 
then initially, the number of the common elements between these two vectors is 
calculated. For example, if Z = 20 and 



Vi =[60 55 52 49 47 43 39 34 33 30 29 22 20 17 14 1 1 9 5 2 l] 



U|=[60 58 55 49 47 41 39 37 33 30 28 25 2017 14 11 9 6 4 2] 
then the number of the common elements is thirteen (13), 

25 Subsequently, it is checked if the number of the common elements between the 

vectors Vj and U| is greater than or equal to the number 0.5 which is called 
"requisite similarity threshold". If, indeed, it is greater than or equal to 0.51*L, we 
proceed to step Ez below. If it is smaller than 0,51*L, then we consider that the set 
of the tests performed so far did not result to a successfiil recognition, so, after 

30 considering U| as the next representative- vector of the model signal, we start the 
comparison procedure again, beginning fi-om the comparison of the vector V| with 
the new Uj . 

El) If the second representative vector of the unknown signal, corresponding 
to the same shift coefficient with V| // , is called and the representative vector 
35 of the model signal corresponding to the sample (^i*/^ ), is called U2, then we 
calculate the number of the conimon elements between these two vectors. 
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Afterwards, we check if the number of the common elements between the 
vectors V2 and U2 is greater than or equal to the "requisite similarity threshold". 
If it is greater or equal, we proceed to step E3 below. If it is smaller, then we 
consider that the set of tests performed so far did not result to a successful 
5 recognition, so, after considering Ui as the next representative-vector of the model 
signal, the comparison procedure starts again beginning from the comparison of the 
vector V| with the new U| . 



10 

E(M.i)) If the (M-1) representative vector of the unknown signal corresponding 
to the same with V| shift coefficient //, is called V^jyj^i) and the representative 

vector of the model signal corresponding to the sample ^^(^^ */ ) is called 

U(M-l) ' ^hen we calculate the number of the common elements between these two 
15 vectors. 

Next, we check if the number of the common elements between the vectors 
%M-l) U(ivi-i) is greater than or equal to the "requisite similarity threshold". 

If it is greater or equal, we proceed to step Em below. If it is smaller, then we 
consider that the set of tests performed so far did not result to a successful 
20 recognition, so, sifter considering Uj as the next representative-vector of the model 
signal, the comparison procedure starts again beginning from the comparison of the 
vector V| with the new . 

Em) If the M representative vector of the unknown signal corresponding to the 
same with V| shift coefficient /, , is called Vjvf and the representative vector of the 

A/-1 

25 model signal corresponding to the sample fi) is called Uivi , then we 

//=1 

calculate the number of the common elements between these two vectors and 
Uivi and we check if it is greater than or equal to the "requisite similarity 
threshold". If it is greater or equal, we proceed to step Em+i below. If it is smaller, 
then we consider that the set of tests performed so far did not result to a successful 

30 recognition, so, after considering as U| the next representative- vector of the model 
signal, the checking procedure starts again beginning from the comparison of the 
vector Vj with the new U j . 

Em^i) First we check how many of the pairs (V|,Ui), (¥2,112),..^ 
(Vm,Ujv|) have, according to the previous comparisons, a number of common 

35 elements in the interval [0.5 1 *L, 0.71 *L]. If the number of these pairs is greater than 
0.34*M, then we consider that the set of tests performed so far did not result to a 
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successfiil recognition, so, after considering Uj as the next representative-vector of 
the model signal, the comparison procedure starts again beginning from the 
comparison of the vector with the new U|. If the number of liiese pairs is 
smaller or equal than 0.34*M then the following check is realised: 

For the pairs of the vectors pairs (V|,Ui), (V2,U2), (Vm^Uj^), 
having already being compared, we calculate the mean value of the number of the 
common elements. If this mean value is greater than or equal to 0.71 *Z. then we 
consider that the comparison between the group of the M representatives of the 
model signal corresponding to the shift coefficient fi that we checked and the 
group of the representatives of the unknown signal, is successful. If the mean value 
is smaller than 0.71 *L then we consider that the set of tests performed so far did not 
resuk to a successful recognition, so, after considering as U| the next 
representative-vector of the model signal, the comparison procedure starts again 
beginning from the comparison of the vector Vj with the new U j 

If all possible vectors of the model signal are unsuccessfully compared with 
one group of representatives of the unknown signal corresponding to the specific 
shift coefficient // , then we repeat the comparison procedure, using the group of 
representatives of the unknown signal corresponding to the next shift coefficient 
yj^l . If the comparison of a specific set of model vectors with all (S+1) groups of 
representatives of the unknown signal is unsuccessful, then we proceed to the 
comparison of the unknown signal with another set of model vectors. 

If the result of the above comparison is successful for a group of the unknown 
signal corresponding to a specific shift coefficient, let*s say we proceed to the 
application of the irrevocable comparison criterion, which will be described below 

As it is already mentioned, the successful application of the first 
aforementioned criterion results to the determination of a group of M 
representatives of the model signal U|» U2»--.9Ujvi which "fit" to the group of 
the representatives of the unknown signal , Vjyj corresponding to the 

specific shift coefficient Since the positions of these vectors in their 

corresponding signals are now known, it is possible to realise a sequence of 
comparisons between vectors of the unknown signal, corresponding to the specific 
shift coefficient with the vectors of the model signal formed at the specific 

positions where the first criterion was satisfied. 

In this way, in the digitised unknown signal of duration from eight (8) to 
sixteen (16) seconds, a window of length is obtained beginning at the unknown 
signal starting point. In this window a fast Fourier transform is applied again and its 
absolute value is obtained. Subsequently, the peaks of the Fourier transform are 
spotted and their positions are multiplied with the shift coefficient which has 
been previously verified that satisfies the first criterion. Then, in each section, the 
peaks are sorted according to their value. In each section to which at least one peak 
has been previously ascribed, the greater peak is kept to form the "representative 
vector of the unknown signal". Next, the L greater value representatives are spotted, 
where the value of L is the same with the one used in the first criterion. The 
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indicators of the sections corresponding to these representatives, sorted in increasing 
order, form a vector that constitutes the "first irrevocable representative-vector of 
the unknown signal". 

Then, the window slides for samples, where the value of is equal to 

5 ^1 ^j^g fluctuates between 30 and 50. For the new window 

D-\ 

position a new vector is calculated, the same way as it was descnbed before, called 
"second irrevocable representative-vector of the unknown signal". The above 
procedure is repeated for over D-2 windows, each one starting at a distance of kj 

, ^, *(A/-1) , 
samples from the start of its previous window, where . ' - 2, 3, 

10 

In this way, finally, a group consisting of D representative- vectors is created 
We will refer to this group with the name "irrevocable group of representatives of 
the unknown signal". 

In order to obtain the final decision if the unknown signal corresponds to the 
15 model signal in hand, the irrevocable group of representatives of the unknown signal 
is compared to elements of the set of the representatives of the model signal, by 
means of a method similar to the first criterion consisting of the steps briefly 
described below: 

TO If the first irrevocable representative-vector of the unknown signal is called 
20 Vj and Uj is called the representative-vector of the model signal corresponding to 
the position, let's say /ii. where the first criterion has been satisfied, then initially we 
calculate the number of the conmion elements between these two vectors. 

Tz) If the second irrevocable representative-vector of the unknown signal is 

called V2 , then this vector is compared with vector \5\ , which is the representative 
25 vector of the model signal corresponding to the position hx-^ k^* f^, where is 
the shift coefficient that has been calculated fi^om the first criterion. 



JO T,D-i)) If the (D-l)"* irrevocable representative-vector of the unknown signal is 

called V(*D_i) and the representative-vector of the model signal corresponding to 

the sample Yu^^j * fs) '^^ ^^^^ ^(D-l)' t^^" calculate the number of the 
y=i 

conmion elements between these two vectors. 

Finally, having calculated the number of the common elements between these 
35 D pairs of vectors, in order to decide for the identification, we check if the two 
conditions stated below are satisfied: 

[Condition 1] At least 0.825 * D fi-om the pairs of the vectors, have common 
number of elements greater than 0.71 * L. 
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[Condition 2] The total number of the common elements of the vectors, 
namely the sum of the common elements of the pairs 

(V;,U; ),(V2 ,U;),...,(Vi-i),UjD-l))X is greater than 0.6875 * D * L. 

If these two conditions are satisfied, then we have successfully recognised that 
the specific musical composition corresponds to the model signal in hand. 

The whole procedure of the identification is described in the Figures 3, 4 and 

5. 
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The method for the automatic recognition of musical compositions and sound 
signals, which is used for the identification of musical compositions and sound 

5 signals played by radio or TV or performed in public places, is based on the 
existence of a procedure, which is applied to the model signals 2md results to the 
extraction of a set of characteristics, which will finally represent each model signal. 
Besides, it is based on a similar procedure, which applies to the unknown musical 
composition or sound signal for the extraction of similar characteristics and, finally, 

10 it is based on a procedure of comparison performed between the representative sets 
of characteristics of the model and the unknown signal. This method is characterised 
by the model sets of characteristics corresponding to division of the frequency 
domain in bands. It is also characterised by two original criteria for the decision of 
the identification, according to which, two musical compositions or sound signals 

15 are identified only when: 

a) A group of M representative vectors of the model signal U| , U2 » Uj^i 

where two successive vectors are calculated at samples having distance , / = 1, 2, 
M-1, "match'' with a group of representatives of the unknown signal Vj , V2, 
Vjyi , which corresponds to a specific shift coefficient . Notice that, the values 
20 of ^y, / = 1, 2, are not necessarily equal, but , in any case, are kept fixed 

throughout the application. The matching between U|, -mUm and V| , V2, 
...» Vj^ is realised by means of the following criterion: 

All comparisons between the vectors of the pairs (V|,U|), (V2,U2), 
( Vj^ , U ) are made and the number of pairs with common elements in the interval 
25 [0.51 * JL, 0.71 * L\ is computed. If it is greater than 0.34 * M then we consider that 
the set of comparisons performed so far did not result to a successful recognition. If 
this number is smaller or equal than 0.34 * M then it is checked if the mean value of 
the number of common elements of vectors of the above pairs (Vj,U|), 

( V2,U2), ( Vjy|,U|vi ) is greater than or equal to 0.71 * L, If it is, then we 
30 consider that the comparison between the group of the M representatives, 
corresponding to the shift coefficient f^y of the model signal in hand and the group 

of representatives of the unknown signal is successful. 

b) A second group of D irrevocable representative-vectors of the model signal 

U|,U2, . .,U^ being calculated at a distance /T/ the one firom its previous, where 
^ * (M - 1) 

35 ki = ^ ^ — ^ , / = 1, 2, ...,Z)-1, which are not necessarily equal, but are, in 

any case, kept fixed throughout the application, "match" with a group of 
representatives of the unknown signal V| , V2, . .,Vd which corresponds to a 
specific shift coefficient , according to the following criterion: 
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• At least 0.825*Z) from the pairs of the vectors 
(V;,U;),(V2,U;),...,(V(*D-i),U(D-l))» have common number of elements 

greater than 0.71 * L. 

• The total number of the common elements of the vectors (namely the one 
5 that results from the summation of the common elements of the pairs 

(V^U;),(V^U;),...,(V(;.,),U*^_,))), is greater than 0.6875 *Z>*Z. 

If both these conditions (a) and (b) are satisfied, then we have successfully 
recognised the specific musical composition. 



0104a70A1 I > 



/ Musical 
/ Composition 



Musical 
Composition 



Digitization 



Digitization 



Windowing 



Windowing 



Feature 
Extraction 
Procedure 



Feature j 
Comparison U- 
Procedure I 



Store Musical 
Characteristics 



Feature 
Extraction 
Procedure 



Results 
Presentation 



End of 
Procedure 



Figure 1 



SUBSTITUTE SHEET (RULE 26) 



wo 01/04870 



Characteristics 
Extraction from 
Model Musical 
Composition 



Choise of the 
I starting Window 



PCT/GROO/00024 



2/5 



Fast Fourier 
Transform 



Adaptive 
Fast Fourier 
Transform 



Values from the 
Discrete Fourier 
Transform 



I Distribution of FFT 
t Peaks to bands 



I Computation of the ; 
I "Representative 
i Vector" of a Window 



Store Vector 



is it the end of the 
Musical Composition 



Musical 
Characteristics 



Noi^ 



New Window 
Selection 



Yes 



/Terminate Procedure for the - 
\specific Musical Composition 



Figure 2 



0104870A1_I_> 



SUBSTITUTE SHEET (RULE 26) 



Choice of the Musical ; 
composition to be i 
compared ! 



I Choice of the first i 
I Window of the 
• Unknown Musical 
Composition 



k=0 



Fast 
Fourier 
Trasform 



i Calculation of the 
j Fourier Transfonn 
i Peaks 



1=0 



! Calculation and Storing 
i of the Representative 
ivector conresponding to f j 



i=i+1 



No 

I 



i>S 



Yes 



Slide Window i 



i k=k+1 



by I, 



J 



I Chose Next 
k>M >-No-^i Window 

i 

— Yes — » . 



Control 



BNSDOCID: <W0 0104870A1J_> 



SUBSTITUTE SHEET (RULE 25) 



wo 01/04870 



PCT/GROO/00024 



'Control 



4/5 



Choise of the first 
stored Musical 
comosttion 



Calculation of the 
position of Vector U. 



{ 1 



i=1 



Calculation of the 
common elements 
between 



V, and U. 



Number of 
Common Elements 
>0,51*L 



Yes 

r 

i=i+1 



l>M 



2 

No ■' / • 



No 



End of Musical 
Comosition ? 



Yes 



Last Stored 

Musical 
Composition 



No 

^ Recognition 

Yes 



No 



No 



No 



No 



Choose next sample of | 
stored Musical I 

Composition i 



Choose New 
stored Musical 
Composition 



Yes 

Number of 
pairs with a few 
common Elements 
>0,34*M 

Yes 



Mean Value of 
Pairs 
> 0.71 *L 



-Yes- 



Final 
Check 



Figure 4 



,0104870A1_L> 



SUBSTITUTE SHEET (RULE 26) 



wo 01/04870 



5/5 



PCT/GROO/00024 



/ Final \ 
\ Chech / 



i=0 



I Choose window j 



Calculate K. 



Fast Fourier 
Transform 



iCalcutationof the representative ; 
' vector that conresponds to a 

specific shift factor \ 



U 

i_ 



i=i+1 



No 



I Calculation of the corresponding 

I representative vector of the model 
I Musical Composition and the 
t number of common elements 



i<D 



^ — Yes^ 



Is condition 1 
true; 



-<~Yes-^ 



No 



Is condition 2 \ ~ 

true; . ' — No^ 2 



Yes 



Succesful 
Recognition 



Figure 5 



SUBSTITUTE SHEET (RULE 26) 

BNSDCXID: <WO 0104870A1_L> 



DTTERNATIONAL SEARCH REPORT 



Internal) i Appllcathm No 

PCT/GR 00/00024 



A. CLASSIFICATION OF SUBJECT MATTER 

"PC 7 GlOHl/12 



According to tntemattonai Patent Classification (IPC) or to both national dassification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system fottowed by classification symbols) 

IPC 7 GIOH GIOL 



Documentation searched other than minimum documentation to the extent that such documents are included In the fields searched 



Bectronic data base consulted during the international search (name of data base and, ¥vhere practical, search terms used) 

EPO-Internal 



C. DOCUMEMTS CONSIDERED TO BE RELEVAWT 



Category * 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to daim No. 



US 5 210 820 A (KENYON STEPHEN C) 

11 May 1993 (1993-05-11) 

column 13, line 29 -column 17, line 58; 

figures 1,2,5,6 



US 5 874 686 A (GHIAS ASIF U ET AL) 
23 February 1999 (1999-02-23) 
column 1, line 57 -column 2, line 3 
line 27 - line 59 
line 31 -column 6, line 59; 



column 2, 
column 4, 
figure 1 



US 5 778 335 A (UBALE ANIL WAMANRAO ET 
AL) 7 July 1998 (1998-07-07) 
column 3, line 10 -column 6, line 2; 
figure 1 



□ 



Further documents are listed in the continuation of box C. 



Patent family members are listed in annex. 



" Special categories of cited documents ; 

"A" document defining the general state of the art which is not 

considered to be of particular relevarx:e 
'E' easier documeru but pubKshed on or after the international 

fifing date 

"L" document which may throw doubts on priority claim(s) or 
which is dted to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P' document put>lished prior to the international fUtng date but 
later than the priority date d aimed 



later document published after the international filing date 
or priority date and not in conflict with the application txit 
cited to understand the principle or theory underlying the 
invention 



document of particuiar relevance: the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 
"Y" document of particular relevance: the claimed invention 

cannot be considered to invdve an inventive step when lt>e 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

document member of the same patent family 



Date of the actual completion of the international search 



27 October 2000 



Date of mailing of the international search report 



06/11/2000 



Name and mailing address ol the ISA 

European Patent Office. P.B. 5918 Patenflaan 2 
NL - 2280 HV Rijswijk 
Tel. (♦31-70)340-2040, Tx. 31 651 eponl. 
Fax: (+31-70) 340-3016 



Authorized offk;er 



Pulluard, R 



Form PCT/)SA/210 (secxxxl sheet) (July 1992) 



: <WO 0104870A1J_> 



INTERNATIONAL SEARCH REPORT 

If., .matlen on patent family members 



Intemar I Appllcatton No 

PCT/GR 00/00024 



Patent document 




Publication 


Patent family 




ruDucanon 


cited in search report 




date 


meniber(s) 




date 


US 5210820 


A 


11-05-1993 


AT 




T 


15-09-1996 








PA 


CUH 1 / J** 


A 


03-11-1991 








DE 


69122017 


D 


17-10-1996 








DE 


69122017 


T 


10-04-1997 








EP 


0480010 


A 


15-04-1992 








ES 


2091328 


T 
1 


ui-1 1-iyyo 








HK 


133697 


A 


24-10-1997 








JP 


5501166 


T 
1 


04-03-1993 








WO 


9117540 


A 


14-11-1991 


US 5874686 


A 


23-02-1999 


NONE 








US 5778335 


A 


07-07-1998 


NONE 









Form PCT/ISA/210 (patent family annex) (July 1 992) 

<WO 0104870A1_I_> 



